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C*~> ! Abstract 
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It has recently been speculated that the solar neutrino mixing angle is connected to 
CO ' the golden ratio ip. Two such proposals have been made, cot $12 = <p and cos 612 = 

ip/2. We compare these Ansatze and discuss a model leading to cos #12 = (p/2 based 
on the dihedral group D\q. This symmetry is a natural candidate because the angle 
in the expression cos 612 = is simply 7r/5, or 36 degrees. This is the exterior angle 
of a decagon and D\q is its rotational symmetry group. We also estimate radiative 
rS I corrections to the golden ratio predictions. 
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1 Introduction 



The question what kind of flavor model underlies the peculiar features of lepton mixing is 
one of the dominating ones in contemporary theoretical neutrino physics. One hopes that 
precision measurements of the flavor parameters will provide hints towards the symmetry 
principle behind the apparent regularities. We will in this paper discuss one intriguing 
example of this line of thought. 

All these issues are linked to the structure of the neutrino (and the charged lepton) mass 
matrix. Typically, the smallness of \U e s\ and the close-to-maximality of 623 are - in the 
charged lepton basis - attributed to the presence of an approximate (i-r symmetry: 




m v = • D E . (1) 



The eigenvector to the eigenvalue D — E indeed is (0, —1, 1) T , but solar neutrino mixing 
is unconstrained by the matrix given above. If in addition to /i-r symmetry the condition 
A + B = D + E holds, then the value sin 2 (9 12 = | is obtained: the infamous tri-bimaximal 
mixing [1] , which dominates the current theoretical literature on lepton flavor model build- 
ing. However, comparing the tri-bimaximal mixing parameters (sin 2 # 12 = |, sin 2 9 13 = 
and sin 2 #23 = |) with the current best-fit, 1, 2 and 3a ranges [2] (very similar results are 
found in [3]) 

sin 2 9 12 sin 2 13 sin 2 #23 
0.304 OiOl O50 

0.288^0.326 < 0.026 0.44-^0.57 (2) 
0.27^0.35 < 0.040 0.39-^0.63 
0.25^0.37 < 0.056 0.36^0.67 

one notes that there is ample room for mixing scenarios other than tri-bimaximal mixing. In 
this short note we consider alternatives to tri-bimaximal mixing and focus on the fascinating 
possibility to link solar neutrino mixing with the golden ratio Lp = Lp 2 — 1 = 1(1 + y/E). 
Two such proposals have recently been made. The first one is [4-6] 

(A): cot9 12 = V sin 2 9 l2 = = —^j= ~ 0.276 . (3) 

1 + ip 5 + V 5 w 

The second possibility is [7] 

(B): cos# 12 = f ^sin 2 # 12 = \(3 - V ) = ~ 0.345 . (4) 

It can be seen that the both predictions lie within the current 2a range0. The possibility 
(A) has first been noted in Ref. [4], and discussed in more detail in [5], where it was also 
mentioned that A§ might be a candidate for the underlying flavor symmetry group. In this 



1 Actually, prediction (A) would lie very slightly outside the 2a range of Ref. [3], which is sin 2 612 = 
0.278 -j- 0.352. 
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spirit a model based on A5, which can lead to cot 612 = ip, has been outlined in Ref. [6]. 
The reason why A 5 is the candidate symmetry is because this group is isomorphic to the 
rotational group of the icosahedron and its geometrical features can be linked to the golden 
ratio. For instance, the 12 vertices of an icosahedron with edge-length 2 have Cartesian 
coordinates (0, ±1, ±y>), (±l,±<p, 0) and (±<£>, 0, ±1). A peculiar feature of cot 9 12 = <p 
is that the angle gives also tan26 , 12 = 2, and this can be obtained from a simple matrix 
proportional to [5] 

m v oc ( \ J . (5) 



x 1 1 t 

This matrix is invariant under a Z2 symmetry generated by [5] 



V5 

where invariance is fulfilled when S T m v S = m v . 



Now consider the second golden ratio prediction cos #12 = ip/2, which corresponds simply to 
9i2 — 7r/5. A mixing scenario based on this value was proposed with a purely phenomeno- 
logical purpose in Ref. [7]. A unified parametrization of both the CKM and the PMNS 
matrix was constructed by choosing in addition to the lepton mixing angle 9 12 = 7r/5 a 
similar expression for the quark sector, namely 9\ 2 = it/ 12. The resulting value for the sine 
of the Cabibbo angle, sm9l 2 = (V3 — l)/y/8, is also a simple algebraic and irrational num- 
ber. The point made in Ref. [7] was that at zeroth order the CKM matrix is a 12-rotation 
with angle 7r/12, while the PMNS matrix is a 12-rotation with angle 7r/5 multiplied with 
an additional maximal (atmospheric) 23-rotation. To correct the 12-angles of the quark 
and lepton sectors to their respective best-fit values, one needs to multiply both zeroth 
order mixing matrices with a small 12-rotation. It turns out that one can achieve this with 
a universal (i.e., the same for quarks and leptons) angle ei 2 — —0.03 [7]. 
In the present letter we concentrate on the possible theoretical origin of the golden ratio 
prediction (B). We stress that flavor models based on the symmetry group D\q are natural 
candidates to generate 6*12 = tt/5. The dihedral group D\q is the rotational symmetry 
group of a decagon and the exterior angle in a decagon is nothing but 7r/5, or 36 degrees. 
Indeed, we will present a model based on D w in the next Section [2J We remark that 
also -D5, the rotational symmetry group of a regular pentagon, could be possible. In a 
pentagon the length of a diagonal is if times the length of a side. The triangle formed by 
the diagonal and two sides has one angle of 108° (the internal angle) and two angles with 
36° each. However, here we focus on D w because it turns out that the vacuum alignment 
we need in our model is simplified due to the larger number of representations in D w . 
Note that just as considering A 5 for the golden ratio prediction (A) was motivated by 
geometrical considerations, the use of the (mathematically simpler) pentagon or decagon 
symmetry group is here motivated by prediction (B). These are examples for the hope 
mentioned in the beginning, namely that precision measurements may give us hints towards 
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the underlying symmetry behind flavor physics^]. 

The present paper is build up as follows: after discussing general symmetry properties of 
mass matrices with cos #12 = cp/2 and an explicit Diq model in Sec. [2] we will in Section [3] 
deal with renormalization group corrections to both golden ratio predictions (A) and (B), 
before we conclude in Section HI 



2 Golden Ratio Prediction 9 12 = n/5 and Dihedral 
Groups 

We have seen in the Introduction that there is a simple Z 2 under which a mass matrix 
generating cot 6*12 = <£> is invariant, see Eqs. (jSJ) and (jSJ). The second golden ratio proposal 
(B) in Eq. (TjJ corresponds to tan26 , i 2 = y/i + (f 2 / (<p — 1), and therefore it diagonalizes a 
less straightforward matrix. Nevertheless, in this case one can make use of Z 2 invariance 
as well, however the charged lepton sector has also to be taken into account. We will first 
discuss this for the simplified 2-flavor case with symmetric mass matrices, before making 
the transition to dihedral groups and then to the explicit model based on D\q that we will 
construct. 



The generators of the Z 2 under which the neutrino mass matrix m v and the charged lepton 
mass matrix have to be invariant are 



S, 



vl 







(7) 



respectively. The matrices m v and are invariant when they have the following structure: 



m, 



B 



vj 



P, 



vl 



e**".< 
1 
A v j B v . 
B u e A u 



Ayl B V J 



B„, A 



1 

e _< *^ 



Qui ■ 



The inner matrix can be written as 

v,t B Vji \ = ^ diag ^ b A Vt i+B v t)U v j, where U vi 

E>vl A v i I 



1 1 

V2 V2 



(9) 



The total diagonalization matrices of m v and mi are \J V ^ = P u> £ U v> i and the physical mixing 
matrix is their product U = U\ U v = U\ P\ P v U v . The 11-element is found to be 



Wel\ 



cos-($ i , - 



(10) 



2 Tri-bimaximal mixing is usually obtained with models based on A4,, the symmetry group of a tetra- 
hedron [8]. Here the angle between two faces (the dihedral angle) is 2#tbm, where sin 2 #tbm = g- 
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The fact that a non-trivial phase matrix lies in between the two maximal rotations U\ 
and U„ is crucial. Obviously, at this stage any mixing angle can be generated. However, 
the observation made in Refs. [9] was that the phase factors in Eq. (0) can be linked to 
group theoretical flavor model building with dihedral groups D n . To make the connection 
from Eq. ( flOl) to dihedral groups, we note that the flavor symmetry D n has 2-dimensional 
representations 2j, with j = 1, . . . , | — 1 (j = 1, . . . , for integer (odd) n, generated by 

e 2 ™« \ , „ / 1 



Z 2 subgroups are generated by 



e - 2m - fc 



BA = I - I (12) 



with integer k. This is just the required form of a Z 2 generator in Eq. (|7j) . It is now possible 
to construct models in which the two fermions transform under the representation 2j of 

_D n , and D n is broken such that m„ is left invariant under B A kv and is left invariant 
under B A kl [9] . Consequently, the relation in Eq. (ITU1) is obtained and we can identify 



Well 



cos 7r — (A;^ — ki 
n 



(13) 



Hence, a natural candidate to implement the requested value of 7r/5 is e.g., D w . This is 
no surprise given the observation that we made in the Introduction, namely that 7r/5 is 
the exterior angle of a decagon and that D\q is its rotational symmetry group. 



We continue with an explicit model: we work in the framework of the MSSM without 
explicitly introducing right-handed neutrinos. Majorana masses for the light neutrinos are 
thus generated by an effective operator coupling to two Higgs vacuum expectation values 
(VEVs). We augment the MSSM by a flavor symmetry x Z^. The symmetry D\q is 
used for our prediction of the solar mixing angle, while the auxiliary Abelian symmetry 
Z 5 separates the charged lepton and neutrino sectors. Due to the flavor symmetry, no 
renormalizable Yukawa couplings are allowed for the charged leptons and the dimension 
5 operator giving mass to the neutrinos also vanishes. Mass for the leptons is generated 
by coupling them to gauge singlet flavons, which acquire VEVs and thereby break the 
flavor group. The charged lepton masses are thus generated by dimension 5 operators, the 
neutrino masses by dimension 6 operators^. 

The transformation properties of the MSSM leptons and Higgs fields, as well as the rep- 
resentations under which the flavons transform, are given in Table [lj The multiplication 
table and the Clebsch-Gordan coefficients of D w are delegated to Appendix [A] Note that 

3 In a model including quarks, this may explain m T <C m t without invoking a large tan/3. 
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Field 


h,2 


h 


e 1.2 




h u ,d 


a e 


Xl,2 


Cl,2 


Pl,2 


a u 




xh 




D w 


2 4 


II 


2 2 


II 


II 


II 




2 3 


2 4 


II 




2 2 


2.3 


z 5 


uo 


uo 


UO 2 


uo 2 


1 


LU 2 


LU 2 


cu 2 


cu 2 


w 3 


to 3 


w 3 





Table 1: Particle content of the D w model: U are the three left-handed lepton doublets, 
ef are the right-handed charged lepton singlets and h Uj d are the MSSM Higgs doublets. 
Furthermore we have flavons a e , xt 2 , £12, Pi 21 ° v ^ ¥121 X12 an d £12 which only transform 
under Di x Z 5 . The phase uo = e~ is the fifth root of unity. 



the fermions and the flavons that couple to them are all in unfaithful representations of 
Dio (i.e., in 2_2 and 2^), so that here a D 5 structure would have sufficed. However, the 
full -Dio structure is needed to achieve the desired vacuum alignment. We can continue by 
constructing the Yukawa superpotential, giving the leading order terms for both charged 
lepton and neutrino masses: 



w Y = y\ (k e c 2 x e 2 + k e\ xl) "r + 2/2 (h e °i Pi + h e c 2 p e 2 ) + y% (h e\ x\ + h e c 2 xt) '' 



A 

h h 
+ yt(he c 3 p e 2 + l 2 e c 3 p e 1 ) — + y e 5 l 3 e c 3 a e j- 



A 



A 
(14) 



h 



h 



h 2 



+ ylhh^^+yVzhcr^+yZihhxl + hhx^Tt+yshk^Ti- 



A 2 



A 2 



A 2 



As we will show below in Appendix [Bj introducing appropriate "driving fields" and mini- 
mizing the flavon superpotential leads to the following VEVs for the flavons: 



(Xl) 

(xl) 



e 5 



<£?> 



3irik 

e 5 



(Pi) 
(P%) 



4-Kik 

e 5 



iirik I , (15) 



where k is an odd integer between 1 and 9, and 



(Xi) 
(X2) 



w,. 



(16) 



The VEVs of the singlet flavons (a e ) = x e and (a u ) = x v are assumed to be also non- 
vanishing. The VEV structure leads to the following mass matrices: 



mi 



(h d ) 
~A~ 



( 



y 2 Z e 



y\ Ve 

2irik 

\y 3 v e e 5 
(K) 2 



2nik 

yfv e e 5 


4-tt ik 

y%z e e 5 


\ 


ylz e e 5 


y%z e 




Vs v e 


yt x e 


J 


v yi 






y2 w u 


■') 

y3%u 1 





(17) 
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To see that indeed the golden ratio prediction is obtained from the above two matrices, 
note that for the choice k = 3 the relevant matrix me ml takes the form 



m e m t 




Be~ 2iq> 
A 

D e -*P+*) 




where $ 



47T 

T 



The quantities A, B, D, G are real and positive, 5 is a phase. To obtain the golden ratio 
prediction for the solar mixing angle, we have to set in Eqs. (I15II17P k = 3 or k = 7. 
From the other possibilities k = 1 or k = 9 would give a solar mixing angle of while 
k = 5 would give a vanishing solar mixing angle. This small number of degeneracies can 
not be resolved by the flavon potential. Looking at the last matrix mem\ in Eq. (jTBI) . 
one immediately recognizes the Z 2 -invariance of the upper left 12-block, which is just the 
invariance we were seeking for, see Eq. (jSJ). To be precise, the D w was broken in a way 
that mg m\ is left invariant under B A 3 , while the neutrino mass matrix m v is left invariant 
under B A = B. Inserting this in Eq. (1131) . where we have to set j = 4 because the 



el 



first and second left-handed lepton doublets transform as 24, we expect \U, 
which is indeed equivalent to an angle of 7r/5. We will explicitly check this in the following. 



COS ^7T 
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Diagonalizing me m\ with the relation Jj\ me m\ lie 
the matrix 



,t 



diag(mg, m 



mi) is achieved with 



Ue = diag(e- ijCI> ,l,e 



-i(3>+<5^ 



/ 



V 










1 





COS 6> 23 

- sin 623 




(19) 



The diagonal phase matrix on the left is crucial. The rotation angle in the 23-axis is given 

9 fo n 

tan 2fl 23 = ■ 7T —^- A , (20) 



G-A-B 1 
and the charged lepton masses are given by 

A - B , m\ T = \[{A + B + G) ±w (A + B - G)] 

where w - 



'1 + 8D 2 /(A + B-G) 2 . 
The neutrino mass matrix is diagonalized via Ulm v U* = mf^, with 



(21) 



( 



\ 









1 / 



p. 



[22) 



The eigenvalues have in general non-trivial phases which are taken into account in the 
diagonal matrix P, and their absolute values are 



mi 



(/0 s 



A 



—\y2Wv - Vl x v \ 



m 2 



A 2 



m Wy + yi x v 



m 3 



A 2 



1 2/3 x h 



(23) 
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We note that the model makes no predictions about the neutrino masses or their order- 
ing. Nevertheless, one can easily convince oneself that the number of free parameters in 
the model is enough to fit the neutrino and charged lepton masses, as well as the large 
atmospheric neutrino mixing angle 623. The model does in general not predict 623 to be 
maximal, which is not an issue given the fact that it is the lepton mixing parameter with 
the largest allowed range. However, maximal mixing is compatible with the model. We 
have 823 = 7r/4 when G = A + B, in which case m 2 = A + Bt V2D and m 2 as in 
Eq. (|2ip . The fact that there is not more predictivity can be traced to the fact that there 
is a comparably large number of flavon fields required in order to make the model work. 
This is the price one unfortunately has to pay if one insists in the rather peculiar value of 
#12. Given the fact that current data allows for this very interesting possibility, one should 
nevertheless pursue the task of constructing models leading to it. 
The final PMNS matrix is 

U = U\ U v . (24) 

One finds that U e s is vanishing and that atmospheric neutrino mixing is governed by 
tan 2# 23 given by Eq. (|20|) . As mentioned above, the PMNS matrix has a non-trivial phase 
matrix including $ in between the two maximal 12-rotations, one of which stems from Ue, 
the other from U„. As discussed above, this is the origin of the required result. Indeed, 
the 12-element of U is 

\U e2 \ 2 = sin 2 $ = sin 2 71-/5, (25) 

and due to U e 3 = this is just sin 2 6*12. We have thus achieved our goal of predicting 
812 = 7r/5. As discussed in Appendix El higher order corrections to the scenario, as well as 
flavor changing neutral currents, can be estimated to give only very small contributions. 



3 Renormalization Corrections to the Golden Ratio 
Predictions 

It is worth discussing renormalization group (RG) effects to the golden ratio predictions, 
because any symmetry leading to the predictions discussed in this paper could presumably 
be operating at a high energy scale A, and the observables have to be evolved down to the 
low energy scale A. Note that RG corrections to \U e3 \ and #23 are typically suppressed with 
respect to the running of 612 by a factor of Amg/Am^. As the initial values of both \U e s\ 
and 823 need not to be specified here (other than being small or close to maximal, respec- 
tively) we do not comment on their RG-shift. We will stay here model-independent and 
estimate the corrections as a function of the unknown neutrino mass values and ordering. 
An expression for 612, where the dot denotes the derivative with respect to t = In fi/fio 
with jj, the renormalization scale, is given e.g., in [10]. One can therefrom estimate the 
shift for the solar neutrino mixing angle: 

0i2^0? 2 + A;i2eBG, (26) 



S 



where 9® 2 is the initial value of 6*12 (here given by Eq. or (J4])) and 

m 2 A 

e R G = c— ^ln- (27) 

with v u = 246 GeV, c = -§ in the SM and (1 + tan 2 0) in the MSSM. Neutrino physics is 
included in 



2ia 



2 



k 12 = sin 2Q\ 2 sin 2 6" 23 \ . (28) 

Consequently from Eq. (|26|) one finds! 

sin 2 9 12 ^ sin 2 6 Q l2 + k 12 e RG sin 2Q\ 2 . (29) 

Note that the Majorana phase a can suppress the running. As well known, 12 decreases 
in the SM and increases in the MSSM, independent on the sign of Am\. The following 
numerical estimates are done with sin 6>2 3 = f , A/A = 10 10 and with Am|, Am 2 , fixed for 
simplicity at their current best-fit values [2]. In the normal hierarchy (NH, ~ a/ Am^, 
m2 — 1/Anig ^> mi) the running in the SM is completely negligible. In case of the 
MSSM, even for tan (3 — 40 the shift in sin 2 # 12 is not more than 1.5 %. This changes in the 
inverted hierarchy (IH, m 2 ~ rrii ~ a/ Am^ 3> 777.3), where in the MSSM and tan/3 = 10 
the value of sin 2 # 12 can increase by around 10 %. In the SM, again, the shift is with less 
than half a percent not measurable. For quasi-degenerate neutrinos with a common mass 
scale of 0.2 eV the SM allows shifts of around 3 %, whereas in the MSSM the shift can be 
as large as the value of sin 2 #12, even for small values of tan/3 = 5. We illustrate this in 
Fig. CE1 where we used Eq. ( 1291) to show the RG-induced shifts of sin 2 #12 for a normal mass 
hierarchy (SM and MSSM with tan/3 = 40), an inverted hierarchy (SM and MSSM with 
tan/? = 10), as well as quasi-degenerate neutrinos (QD, smallest mass 0.2 eV for the SM 
and MSSM with tan/3 = 5). In case of a normal and inverted hierarchy we have chosen (at 
high scale) 0.001 eV for the smallest neutrino mass. To a good approximation and unless 
in the MSSM tan /3 is very large, the running of the neutrino masses can be described by a 
rescaling, with basically no dependence on the other neutrino parameters [10]. Because k\2 
from Eq. (1281) has the masses appearing in the denominator and numerator, their running 
cancels in our approximation as long as IA^crcI 1- The range of the corrections in 
Fig. [T]is due to the unknown Majorana phases. For illustration, we also include the shifts 
for tri-bimaximal mixing. 

To bring 6 I 12 very close to the best-fit value, the prediction (A) requires the MSSM and IH 
or QD, while prediction (B) (and tri-bimaximal mixing) requires the SM with rather large 
neutrino masses. If future data leads to more precise determinations of sin 2 Q\ 2 and other 
neutrino parameters, one will be able to rule out some of the existing possibilities. 



4 Inserting sin 2 9® 2 — 5 and sin 2 #23 = \ in the following and the last expression reproduces the results 
fromRef. [11]. 
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QD+MSSM 
QD+SM 

IH+MSSM 
IH+SM 

NH+MSSM 
NH+SM 
initial 
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I | I I I ] | I I I I | I 



I | I I I | I I I I | 



i i i i I i I i i I i i i 



I I I | I I 



I | I I I I | I I I I | I I 



• 


best-fit 




■ 


TBM 




▲ 


cot 9 [2 = 


<P 


T 


cos B 12 = 


<p/2 




2a 






la 





I I I I I I II I I I I I II I I I I I I I I I I I I I I 



0,27 0,28 0,29 0,3 0,31 0,32 0,33 0,34 0,35 0,36 0,37 



sin 6 



12 



Figure 1: RG-induced shifts on sin 9\2, estimated from Eq. ( 1291) . for the two golden ratio 
proposals and for tri-bimaximal mixing (TBM). The current best-fit value, as well as the 
1 and 2a ranges are also indicated. In case of the MSSM we have taken tan/? = 40 for 
a normal hierarchy (NH), tan/3 = 10 for the inverted hierarchy (IH) and tan/3 = 5 for 
quasi-degenerate neutrinos (QD) with a mass scale of 0.2 eV. The line for the SM and a 
normal hierarchy cannot be seen because the effect is too small. 

4 Summary 

Precision flavor data may give hints towards the underlying physics. We have stressed in 
this paper that current data implies that the golden ratio ip can be connected to solar 
neutrino mixing. With cot 6*12 = <p and cos #12 = (p/2 there are two appealing possibilities, 
not too far away from current best-fit values and compatible with current 2a ranges. We 
have compared these values, estimated radiative corrections and in particular proposed a 
model based on the dihedral group Diq leading to the relation cos 612 = <p/2. The angle 
leading to cos #12 = (p/2 is 9\ 2 = 7r/5 and closely linked to the symmetry of a decagon, 
which naturally leads one to consider its rotational symmetry group Diq. 

Acknowledgments 

This work was supported by the ERC under the Starting Grant MANITOP and by the 
Deutsche Forschungsgemeinschaft in the Transregio 27 "Neutrinos and beyond - weakly 
interacting particles in physics, astrophysics and cosmology" (W.R.). A.B. acknowledges 
support from the Studienstiftung des Deutschen Volkes. 



10 



X 


ll 


±2 


13 


14 


ll 


ll 


12 


13 


14 


12 


12 


ll 


U 


13 


13 


13 


u 


ll 


12 


U 


U 


13 


12 


ll 



X 


2.1 2 2 


2.3 


2 4 


2.1 


ll + 12 + 2 2 ll + 2 3 


2-2 + 2.4 


I3+I4 + 23 


2 2 


2i + 2 3 li+l 2 + 2 4 


13+14 + 21 


2-2 + 2.4 


2 3 


2-2 + 2.4 I3 + I4 + 2_i 


ll+l2 + 2 4 


2l + 2 3 


2 4 


I3 + I4 + 23 2 2 + 2.4 


2.1 + 2.3 


ll + 12 + 2 2 



ll,2x2j = 2 j; l 3 ,4x2 j = 2 



Table Al: Multiplication rules for the dihedral group D w , which has four two-dimensional 
and four 1-dimensional irreducible representations. 



A Multiplication Rules and Clebsch-Gordan Coeffi- 
cients of Diq 

We present here the Clebsch-Gordan coefficients for D^. The multiplication rules for the 
Kronecker products are given in Table lATl For s, ~ lj and (ai, a 2 ) T ~ 2j we find 

I 31 " 1 )~2: , ( S2fll )~2* , ( S3a2 )-2 5 , and 
^ Sl a 2 J -J V -s 2 a 2 J -J \ s 3 ai / ~° J 

The Clebsch-Gordan coefficients for the product of (ai,a 2 ) T with (bi,b 2 ) T , both in ~ 2j, 
read 

a±b 2 + a 2 bi ~ 1^ , a x b 2 — a 2 b\ ~ I2 , 

( aih ]~ 2- or ( a2h U2. 

y a 2 & 2 / _J V ai&1 / _J ' 

depending on whether i = 1, 2 or i = 3, 4. For the two doublets (ai,ci2) T ~ 2j and 
{bi, b 2 ) T ~ 2: we find for i + j 7^ 5 
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Field 




< 2 


tOe 
Sl,2 




x?r 2 


SI, 2 


z 5 


13 

to 




2 3 

a; 


14 

u 4 


cu 4 


2 3 

cu 4 



Table Bl: Transformation properties of the driving fields under D w x Z 5 . Again u> 
is the fifth root of unity e~ . 



If i + j = 5 holds the covariants read 

aib\ + a 2 b 2 ~ I3 , ai&i - ^2^2 ~ I4 , 
( aib 2 \ ^ 2 or ( fl2&1 ^ ~ 2 

V a 2b\ J _k V a l & 2 / _k ' 

Again, the first case is relevant for k = i — j, while the second one is valid for k = j — i. 

B VEV Alignment of the D w x Z 5 Model 

To obtain the necessary vacuum alignment in the flavon potential, we need to introduce 
a U(1)r and driving fields [12]. Regular .R-parity is a subgroup of the U(1)r. To ensure 
a supersymmetric Lagrangian, the superpotential must have a U(1)r charge of 2. The 
superfields containing the SM fermions have an i?-charge of 1, while the Higgs fields have 
an i?-charge of zero. Hence, for the Yukawa superpotential given in Eq. (fl4"l) to be viable, 
the flavons also need to have a vanishing i?-charge. Consequently, for the flavon super- 
potential one needs to introduce additional flavor- charged fields, having an i?-charge of 2. 
The transformation properties of these driving fields are given in Table IB 1 1 The flavon 
superpotential can then be divided into two parts 

W f = W f,e + Wf,u , (Bl) 

where Wf^ e and Wf jV are responsible for the vacuum alignment of the flavons contributing to 
the charged lepton and neutrino masses, respectively. We begin by considering the charged 
lepton part: 

Wf,e = a e (Xt Cl + X2 £2) ^ + b e (Xt £2 Vl + X2 £l + Ce (Cl P% + tl P\ <P?) 

+ d e (C2 sT + Ci e/) ° e + fe (Zl p\ £i° e + C 2 P e 2 ■ (B2) 

As the flavor symmetry is broken at a high scale, the scalar potential can be minimized in 
the supersymmetric limit. The flavons and driving fields are not charged under any gauge 
group, so the scalar potential is given by the F-terms alone. Hence, we can determine the 
supersymmetric minimum of the potential by setting the F-terms of the driving fields to 
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= 0, 


be Xt C2 + CeClp e 2 


= 0, 


K XlCl + CeQpt 


= 0, 




= 0, 


deClCT e + feC2P2 


= 0. 



zero: 

dwf te 

Q^Oe 
dWf te 

d^ e 

dWf te 

~w 

Similarly, from the neutrino part 

Wf,» = ^(Xiei-^^)^ + ^(^Xi^ + ^Xa^) + C,(e^ + ff^)^ (B3) 
+ d u (tf C2 Xi + V2 Ci Xt) + fu mf xt + Xx) + 9u (X2 XT + Xi Xt) ° v , 

we obtain a minimum of the potential by setting the F-terms of the driving fields to zero: 

dxt 

® W f> v A ,„v tv 1 t ( ,Jf\2 



dxt 
dw ftU 



d^Ui+fArt) 2 + g»x> u = o 
Krtxi + CvCiv u = 0- 



7 = u »f 

As advocated above, these two sets of equations are uniquely solved by the VEV configura- 
tions given in Eqs. (jT5l) and (fISl) . where we have set a possible relative phase in the doublet 
of VEVs of the flavons in the charged lepton sector to zero. This can be done without loss 
of generality, as only the phase difference between the two sectors is phenomenologically 
relevant. We have also assumed that none of the parameters in the superpotential vanish. 
For the charged lepton sector, the flavon VEVs w e and x e are free parameters (which we 
take to be non-zero), while 

4nik C e (I e X e 8-nik d e X e , * 

Similarly v v and x v are free parameters (again taken to be non-vanishing) and 

w - CvfuWl _ Kf v vl 

Cu 9v xi - K d v vf, c v g v x l v - b v d v v* 
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The driving fields themselves are only allowed vanishing VEVs, as can be inferred from 
considering the F-terms of the flavons. Note, that since we can not make the cutoff scale 
A arbitrarily large, we need to take into account NLO corrections to both the Yukawa and 
flavon superpotentials. We also should be careful in what regards potentially dangerous 
flavor changing neutral currents induced by the flavons. All this could be taken into account 
by carefully studying the mass spectrum of the scalars. Given the sizable number of fields 
this is a formidable task, but fortunately it suffices to make some general estimates, which 
agree well quantitatively with a lengthy explicit calculation in a similar model [13]: the r 
lepton mass, see Eq. ( l2Ti) . is of order (/) v/A, where (/) is a flavon vev, v the Higgs vev 
(~ 10 2 GeV) and A the cutoff scale. The neutrino mass, see Eq. ( 123|) . is of order (/) v 2 /A 2 . 
With the charged lepton r mass ~ GeV and the neutrino mass ~ 0.1 eV it follows A ~ 10 12 
GeV and (/) ~ 10 10 GeV. Now we can estimate that the flavon mass is also of order of 
(/). NLO corrections to the potential, and therefore to the neutrino and charged lepton 
mass matrices, are of order (f)/A ~ 1CT 2 and therefore under control. Any potentially 
dangerous flavor changing neutral currents are also suppressed by the heavy mass scale 

(/)■ 
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